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Levy flights constitute a broad class of random walks that occur in many fields of research, from 
animal foraging in biology, to economy to geophysics. The recent advent of Levy glasses allows 
to study Levy flights in controlled way using light waves. This raises several questions about the 
influence of superdiffusion on optical interference effects like weak and strong localization. Super 
diffusive structures have the extraordinary property that all points are connected via direct jumps, 
meaning that finite-size effects become an essential part of the physical problem. Here we report on 
the experimental observation of weak localization in Levy glasses and compare results with recently 
developed optical transport theory in the superdiffusive regime. Experimental results are in good 
agreement with theory and allow to unveil how light propagates inside a finite-size superdiffusive 
system. 
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Levy flights are maybe the most general class of ran- 
dom walks, of which the commonly known Brownian mo- 
tion is a limiting case Q • They are governed by Levy 
statistics 0, which have the fascinating property that, 
depending on the value of one control parameter, they 
can exhibit a diverging variance [^j . This leads to a phe- 
nomenon called superdiffusion, being a diffusive process 
in which the mean square displacement increases faster 
than linear in time [a, Levy flights are common in 
nature and appear, for instance, in animal food searches 
0, II] , laser cooling of cold atoms 9], evolution of the 
stock market [lo| , astronomy JT,] , random lasers ^ and 
turbulent flow 131. 



The recent development of Levy glasses [ij] and care- 
fully prepared hot atomic vapours 'l5|, have allowed the 
observation of Levy flights of light waves and the result- 
ing superdiffusion process. Since interference effects play 
a dominant role in light transport, this raises the natural 
question how interference influences optical superdiffu- 
sion - a concept which has not been addressed so far. In 
regular - diffusive - disordered optical systems, interfer- 
ence leads to speckle correlations, and weak and strong 
localization effects, which all have been studied exten- 
sively over the last two decades [3 121- Localization, 



in particular, leads to a complete halt of transport that 
can confine light waves in random patterns. Since the su- 
perdiffusion induced by Levy statistics tends to enhance 
transport, one would expect it to counter- act localization 
induced confinement. 

Among all interference phenomena in random optical 
materials, maybe the most robust is that of weak lo- 
calization J^]. It is observed in the form of a cone of 
enhanced backscattering, which contains information on 
the path length distribution deep inside the random sys- 



tem, and which has been observed in recent years from 
several diffusive random structures Some of us 

have recently shown theoretically that weak localization, 
or coherent backscattering, can be observed from Levy 
glasses and that a superdiffusive approximation can be 
used to predict its behaviour [26]. In this paper we re- 
port on the experimental observation of weak localization 
from superdiffusive materials, which constitutes the first 
observation of an interference effect in transport based on 
Levy statistics. We find a good agreement with superdif- 
fusive transport theory and show how the backscattering 
cone can be used to extract the Green's function in a 
Levy glass. Contrary to regular diffusive media, in a Levy 
glass, light from an arbitrary depth inside the medium has 
a nonvanishing probability to couple directly to the sur- 
rounding environment. This latter property makes light 
scattering from Levy glasses complex, and has important 
consequences for its (back) scattering properties. 

The Levy glasses under investigation are made of 
jammed-packed microscopic glass spheres, whose diame- 
ter (0) varies almost over two order of magnitudes (from 
5 to 230 pm) following a power-law distribution p(0) oc 
0-(^+i)^ with (3 adjustable parameter. These spheres are 
embedded in a polymeric matrix which matches their re- 
fractive index (n = 1.52) and in which Ti02 nanopar- 
ticles (average diameter 280 nm) have been dispersed 
(see Fig. 1(a)) jl^l- Because of their refractive index 
(n = 2.4) higher than the polymer, these nanoparticles 
act as point scatterers which are not homogenously dis- 
tributed throughout the sample due to the presence of the 
glass spheres (see Fig. 1(b)). As a result, light transport 
is dominated by the long "jumps" that light performs 
propagating through the microscopic spheres. The step 
length distribution that light performs in Levy glasses fol- 
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FIG. 1: (a) Electron micrograph of the interior of a Levy 
glass, (b) Sketch representing the optical mechanism lay- 
ing behind the coherent backscattering cone in a Levy glass. 
I/=300 pm is the thickness of the sample and the diameter 
of the spheres. 



lows a power-law decay as p{l) oc where a is re- 

lated to /3 as a = /3 — 1 for an exponential sampling of the 
diameter distribution 27[. By controlling the diameter 
distribution of the spheres in a Levy glass we can control 
a and, thus, the degree of superdifFusion of the material. 
For a > 2 the system is diffusive, whereas for < a < 2 
the system is superdiffusive. However, in real systems 
the finite size of the sample truncates the step-length 
distribution to the largest sphere diameter, enriching the 
physics behind the transport properties 28-33]. In pre- 
vious publications 14;., 21] Levy glasses were fabricated 
between two microscope slides. In contrast, in this work 
we remove one of the two slides to reduce undesirable 
reflections which affect the quality of the measurements. 
Moreover, the thickness of the sample is approximately 
70 pm more than the largest sphere, which means that 
the step length distribution is truncated on a length scale 
that is slightly below the sample thickness. 

The setup employed follows a standard scheme for co- 
herent backscattering experiments. Light emitted by a 
HeNe laser (@632 nm) is expanded to coUimated beam 
of 1 cm in diameter to ensure a high angular resolu- 
tion of the system. A beamsplitter is used to separate 
the backscattered light from light impinging on the sam- 
ple. Subsequently, the backscattering cone is imaged on 
a CCD camera and the use of a polarizer ensures that 
we observe only the polarization conserving channel [l7| . 
The angular resolution of the setup has been optimized 
to properly investigate the different shapes of the mea- 
sured cones. During the data acquisition the sample is 
nutated to average over different realizations of disorder. 

In order to verify the quality of the setup we fabricate a 
set of diffusive samples made of Ti02 and polymer (with- 
out glass spheres) and measured its backscattering cone. 
The amount of polymer was adjusted to have the same 
average (over the total volume) density of scatterers in 
both diffusive and superdiffusive samples. A cross-cut of 
the measured diffusive cones is shown in Fig. 2(a) to- 



gether with a fit (red curve) obtained from the diffusion 
theory [iH. The retrieved transport mean free path is 
i* ~ 19 pm. In the inset the residuals are displayed as 
blue dots. The good match between experimental data 
and fit shows the quality of the experimental setup. 

We measured the backscattering cone on two different 
sets of superdiffusive samples characterized by a = 1.5 
and a = I. The experimental results are shown in Fig. 
2(b) and (c). By comparing these figures one can notice 
the rising of the tail of the cone as the degree of superdif- 
fusion increases, i.e. when a decreases, following a trend 
which has been theoretically predicted in Ref. [2^. To 
ensure that the Levy cones are not merely broader diffu- 
sive cones we made a best fit using diffusion theory. The 
result is shown in Fig. 2(b) and (c) as red curves and 
the residuals of such fits are shown in the inset of Fig. 
2(a). It is clear that the diffusion theory cannot properly 
describe the optical properties of Levy glasses. 

The coherent backscattering cone for the superdiffusive 
samples under investigation can be calculated by taking 
advantage of the fractional derivative approach developed 



in Ref. [26|. The transport of light in a finite, translation- 



ally invariant, superdiffusive system for a point source 
at xq is described by the stationary fractional diffusion 
equation [26l |: 



Da (V^ - kl) f (x, xo, ki) = -Six - xq), 



(1) 



where V" is the symmetric Riesz fractional derivative 
with respect to spatial derivatives, is the in-plane 
component of the wavevector in free space, /(x,xo,k^) 
is the intensity propagator and Da is a generalized dif- 
fusion constant. The spatial nonlocality of V" makes 
the definition of boundary conditions non trivial [3l| . In 
order to model physical systems, such as Levy glasses, 
the fractional Laplacian operator can be represented by 
an M X M matrix, whose eigenvalues Xi, rescaled as 
Ai — i> Xi{M/L)°' with L the slab thickness, and eigen- 
vectors ipi converge to those of the continuum operator 
as M goes to infinity ^3^. Absorbing boundary condi- 
tions can then be implemented by reducing the infinite 
size matrix to a finite size matrix and Eq. [1] be solved 
by eigenfunction expansion. The knowledge of the in- 
tensity propagator of the system then makes it possible 
to calculate interference effects in a "superdiffusion" ap- 
proximation. In particular, considering a planewave at 
normal incidence on the slab interface and in the Fraun- 
hofer regime, the coherent component of the albedo is 
given by the following expression: 



Aaie) oc - ^ F{xi,X2,t 



'ipi{xi)'ipi{x2) 



(2) 



wherei^(xi,a;2,6') = P{xi)P{x2)P{xi/ cose)P{x2/ cos0) 
describes the attenuation for the amplitude of the inci- 
dent and emergent planewaves in the scattering medium. 
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FIG. 2: (a), (b) and (c) Measured coherent backscattering cone from diffusive samples and Levy glasses with a = 1.5 and 
a = 1, respectively. In red the fit to the experimental data according to the standard diffusion theory. In the insets, residuals 
of the fits. 



9 is the angle between the incident and emergent 
planewaves (see Fig. 1(b)) and k± — |kj^| ~ (27r/A) 0, 
at small angle 9. The amplitude attenuation P{1) was 
modelled as a Pareto-like distribution P{1) = 1 for 
< ; < and P{1) = {Ic/l)^"^^^^^ for / > h, where h 
is the cut-ofF length, as to closely follow the step length 
distribution of real Levy glasses l28| . Internal reflections 
were neglected. 

The results are shown in Fig. 3 for a — 1.5 and a — 1, 
where the only adjustable parameter used is Ic- The in- 
set shows the residuals between theory and experiment, 
which are greatly reduced with respect to the diffusive fit 
(Fig. 2(b) and (c)). Due to the very long tail of the cones 
for Levy glasses and the lack of an analytical expression 
for it, the experimental incoherent background was set to 
the one obtained semi-analytically. It must be pointed 
out that for the calculations we employed a step-length 
distribution which was not truncated, in contrast with 
the real system. This is due to the fact that the propa- 
gator / {x,xo,'k±) of the system has been calculated by 
considering the sample sa translational invariant in the 
in-plane direction and finite in the longitudinal direction. 
The very good agreement between experiment and cal- 
culation shows that the fractional diffusion approach can 
properly describe light interference effects due to multiple 
scattering in the superdiffusion approximation. 

The gentle slope of the cone tails arising from a Levy 
glass is a hallmark of anomalous transport of light. This 
is a striking effect since generally the width of the cone 
is expected to be inversely proportional to the average 
step length £, whereas in a Levy glass the average step 
length £a increases when a decreases. The knowledge of 
the peculiar Green's function inside this medium, which 
dictates the shape of the cone, might shed light on this 
phenomenon. The behavior of such a propagator as a 
function of the degree of superdiffusion can be retrieved 
by a Fourier analysis of the experimentally observed co- 
herent backscattering. In diffusive media the coherent 
albedo can be approximated with the propagator of the 
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FIG. 3: Comparison between the calculated superdiffusive 
cones obtained with the fractional derivative approach and 
the measured Levy cones. 



system in reciprocal space. The propagator is taken for 
a point source at a distance of a mean free path away 
from the boundary {Ac oc f{x = £,xo = £,k±)) [s^. We 
expect a similar Fourier analysis on the coherent albedo 
from Levy glasses to provide information on the intensity 
distribution inside a superdiffusive media. In Fig. 4(a) 
the normalized Fourier transform of the measured and 
calculated coherent albedo as a function of the in-plane 
displacement p = \{'ri — 're)± \ are shown. These curves are 
found to be in good agreement, in particular at relatively 
small values of p, and show a remarkable reshaping as a 
function of a. In Fig 4(b) and (c) the normalized Fourier 
transform of the theoretical cone and the normalized in- 
tensity distribution f{x — la,XQ — £a, p) calculated from 
the fractional diffusion approach, respectively, are shown 
for long p. The qualitative agreement between these two 
figures is evident, in particular in their dependency on 
the degree of superdiffusion. The spatial distribution of 
the propagator is dictated by the power-law step length 
distribution in Levy glasses, which allows light to cou- 
ple directly to the surrounding environment from any 
depth inside the sample. This spatial non-locality ap- 
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FIG. 4: (a) Amplitude of the Fourier transform of the measured and calculated coherent backscattering (CB) for small p. (b) 
and (c) Amplitude of the Fourier transform of the calculated cone and calculated intensity distribution f{x — £a,XQ = £a,p) 
for large p. 



plied to an open system such as a Levy glass results in 
a strong modification of the shape of the propagator as 
a function of a [2^. Firstly, the more cusped shape of 
the f{x = ia,xo = ia,p) leads to a smaller in-plane dis- 
placement and thus to higher tails in the coherent albedo. 
Secondly, the gentle tail of the propagator, which induces 
a cusped cone [2^ , is a consequence of the fact that light 
can escape more easily from its local environment to- 
wards large distances. 

In conclusion, we have reported on the experimental 
observation of weak localization from Levy glasses and 
found strong deviations from the prediction of the dif- 
fusion theory. The recently developed semi- analytical 
model for superdiffusive system is able to reproduce 
the experimental results with a high degree of accuracy. 
From the study of the coherent backscattering cone we 
can retrieve the behavior of the Green's function inside 
real systems in which the long steps probability is non- 
vanishing, not only in Levy glasses but possibly also in 
foams, clouds or strongly heterogeneous porous media. 
Further coherent backscattering studies oriented to in- 
vestigate the influence of real systems parameters (e.g., 
finite size, truncation, quenching) on the superdiffusive 
backscattering and their fingerprints in the dynamics of 
the backscattering cone can be performed to unravel the 
optical properties of real superdiffusive systems. 
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